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Abstract. We review the techniques of lattice QCD calculations for excited hadrons with light quarks and 
outline the future challenges that are faced in calculations with fully dynamical fermions. 

PACS. 11.15.Ha Lattice QCD 



1 Introduction 

In the last 20 years the lattice approach has evolved into a 
powerful quantitative tool for obtaining non-perturbative 
results from ab-initio QCD calculations. A revolutionary 
breakthrough was the understanding of chiral symmetry 
on the lattice, which may be implemented by using a lat- 
tice Dirac operator that obeys the Ginsparg- Wilson equa- 
tion [1]. Thus lattice QCD is now conceptually ready for 
calculations with fully dynamical light quarks, and as al- 
gorithms and computer technology are evolving rapidly, 
more challenging problems are considered. 

An prominent example for this type of more advanced 
calculations is the analysis of excitations of hadrons. Ex- 
cited states are much harder since they, as we will outline 
in the next section, appear only in sub-leading terms of 
the Euclidean correlators one analyzes on the lattice. Al- 
though powerful methods are available for tackling this 
problem, it is still a challenge to construct hadron inter- 
polators that have a strong overlap with the excitations, 
and at the same time may be implemented numerically in 
a cost-efficient way. 

About 5 years ago several groups started to systemat- 
ically explore hadronic excitations using lattice QCD [2]. 
A strong motivation for this enterprise is, e.g., the un- 
resolved nature of the first positive parity excitation of 
the nucleon, the N(1440) Roper state. Although calcula- 
tions at relatively large volumes with light quarks were 
performed, there is still no definitive answer to the Roper 
puzzle from the lattice. In particular the role of dynamical 
sea quarks still has to be understood more clearly. Nev- 
ertheless, the Roper problem has considerably pushed the 
development of the techniques for excited hadrons. 

After a brief introduction to lattice methods, in this 
contribution we review the recent developments and out- 
line how a state of the art calculation for excited hadrons 



with light quarks proceeds. We illustrate the current sta- 
tus by discussing a few selected results in more detail. 

The above mentioned progress towards fully dynami- 
cal simulations with light quarks gives rise to new chal- 
lenges: So far most of the calculations were done in the 
quenched approximation where excited states cannot de- 
cay. The same is true for a dynamical calculation when the 
quarks are heavy such that decays are not possible due to 
kinematical reasons. However, slowly data from simula- 
tions with light dynamical quarks become available and 
the problem of particle decay and scattering states has to 
be faced. In the continuum a scattering state has a contin- 
uous spectrum of energies and for a particular value of the 
relative momentum its energy may be degenerate with the 
energy of a bound state. In the Euclidean correlator one 
studies on the lattice, this degeneracy leads to a mixing of 
scattering and bound states. The tools for disentangling 
the two types of states are ready in principle - one uses a 
finite volume analysis - but the implementation of these 
ideas in full lattice QCD calculations is still in its infancy. 
We briefly address these problems, together with other fu- 
ture challenges, such as the evaluation of matrix elements 
for excited states. 



2 Excited states in lattice QCD 

2.1 Euclidean correlators on the lattice 

Before we come to discussing excited states on the lat- 
tice, we need to review the central tool in lattice QCD, 
Euclidean correlators. In any lattice calculation one eval- 
uates numerically the correlators of some operators O. The 
Euclidean correlator of two operators Oi, O2 is defined as 
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In this expression H is the QCD Hamihonian and t and 
T are real numbers, often referred to as Euclidean time. 

The interpretation of the Euclidean correlators is ob- 
tained with the help of the (unknown) eigenstates \n) of 
the Hamiltonian which obey H\n) — En\n), where En de- 
notes the energy of the eigenstate \n). If one uses them to 
compute the traces in ([T]) and ([2]), and inserts the unit in 
the form 1 = J^n \^){M ^^ the left of Oi in ([1]), one finds 
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In this spectral representation for the Euclidean 2-point 
function, the sum runs over all eigenstates \n) of the QCD 
Hamiltonian. Each term comes with a Boltzmann factor 
containing the corresponding energy En- The energies are 
normalized relative to the vacuum state |0) with energy 
Eq = 0. The Boltzmann factors are multiplied with matrix 
elements of the operators Oi between the vacuum state |0) 
and a physical state \n). 

The spectral representation ([2]) is a powerful tool. If 
one, e.g., wants to compute the mass of the nucleon, one 
uses for both Oi and O2 an operator On with the quantum 
numbers of the nucleon. For this choice the matrix element 
(njOJYlO) is non- vanishing only for those terms in the sum, 
where \n) is the state \N) of the nucleon or one of its 
excitations |iV'), |iV") ... . Thus we obtain: 
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Obviously we can extract the energy of the nucleon from 
the exponential decay of the Euclidean correlator. Fur- 
thermore, the overall factor is related to the matrix ele- 
ment of our operator with the physical nucleon state \N). 
Having convinced ourselves, that the Euclidean cor- 
relator can be used to extract energies and physical ma- 
trix elements, we need to discuss the way the Euclidean 
correlators are evaluated in lattice QCD. Instead of the 
unknown eigenstates of the QCD Hamiltonian, one uses 
eigenstates of the field operators for computing the trace 
in ([T]). The exponential function of the Hamiltonian is 
treated with the help of the Trotter formula, which in- 
troduces an auxiliary direction that in the end gives rise 
to a 4-dimensional lattice formulation for the Euclidean 
correlation function, given by 
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The Euclidean correlation function is represented as a 
path integral over all degrees of freedom, the gluons G 
and the quark fields q, q. All fields live on a 4-dimensional 
lattice, and the path integration 2?[G', 5, g] is implemented 
as the product of the integrals over the classical degrees 
of freedom on all lattice points. Each configuration of 
the fields G, q, q is weighted with the Boltzmann factor 



exp(— S*), where S[G,q,q\ is a lattice discretization of the 
QCD action (derivatives are replaced by finite differences 
on the lattice etc). In the path integral the operators Oi 
appear as monomials Oi of the classical field variables 
G,q,q. We will refer to these monomials as interpolators 
from now on. The Euclidean time arguments t and {t =) 
determine from which time-slice of the 4-dimensional lat- 
tice the field variables are taken for building the interpo- 
lators. 

The expression ([S|) has the form of an expectation 
value for a statistical system in the canonical ensemble. 
Thus the numerical methods from statistical mechanics, 
in particular Monte Carlo simulation techniques, can be 
taken over to lattice QCD. Using these we can evaluate 
the Euclidean correlators numerically and extract the re- 
sults for energies and matrix elements using the spectral 
representation ^. However, we stress at this point, that 
the Monte-Carlo methods which are used for evaluating 
the Euclidean correlators give rise to statistical errors. As 
we will see below, this limited accuracy implies that spe- 
cial methods for extracting the physical observables from 
the correlators have to be developed. 



2.2 Why are excited states so difFicult? 

Having outlined the basic steps of a lattice calculation, 
we can return to the problem of excited states. We have 
discussed that by selecting suitable interpolators with the 
quantum numbers J, J, P, C we can project onto the states 
we want to analyze, and the example of Eq. ^ illustrates 
how to study the positive parity nucleon channel. How- 
ever, when we want to consider the first positive parity 
excitation of the nucleon, the N(1440) Roper state, we 

find that this state has the same quantum numbers ^ ^^ 
the nucleon ground state. Thus all we can obtain from our 
Euclidean correlator (ONit)Olf{0)} is a sum of exponen- 
tials 
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where |A^') denotes the first excited state (with positive 
parity) and M^' is the corresponding mass (from now on 
we use interpolators projected to vanishing momentum, 
such that the energy reduces to the mass, Ei = Mi). The 
dots indicate the contributions of higher excitations with 
positive parity, such as the N(1710). We stress that nega- 
tive parity states, e.g., N(1535) and N(1650), do not con- 
tribute in the tower of excitations, since the parity of our 
lattice interpolator On was chosen positive. 

Eq. ([7|) illustrates clearly why excitations are much 
harder to extract: The spectral decomposition of the Eu- 
clidean 2-point function is dominated by the exponential 
decay from the lightest mass Mn- Only the sub-leading 
term contains the mass Mn' of the first excitation we want 
to study and this contribution is suppressed exponentially 
for increasing t with a factor of exp(— (Ma?/ —Mn) t). This 
implies that only for small values of t we can hope to see 
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a signal of the first excitation. For higher excitations the 
situation is even worse. 

Ahhough the signals of the excited states are sup- 
pressed exponentially, several different techniques for ex- 
tracting their spectrum have been applied to the Euclidean 
correlators. In the next section we will discuss in detail 
the variational method [ro", 11 . However, also other ap- 
proaches have been tried [12 . These attempts employ ad- 
vanced fitting techniques or other analysis tools for getting 
both the leading and the sub-leading exponential or even 
try to reconstruct the spectral density of the two point 
function. The results typically suffer from the statistical 
errors which make multi-exponential fits rather unstable 
and from the fact that excitations produce significant con- 
tributions only for small t. Consequently only a limited 
amount of information is available. In general the impres- 
sion is that the alternative methods [12] have not reached 
the quality of results from the variational method which 
we discuss in detail in the next section. 
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Fig. 1. Euclidean correlators for different interpolators with 
the quantum numbers of the pion. We plot the correlators as a 
function of Euclidean time t, using a logarithmic scale for the 
vertical axis. 



3 Lattice technology for excited states 

3.1 The variational method 

In the last section we have addressed the problem of ex- 
tracting excited state masses from the sub-leading expo- 
nentials of Euclidean correlators. The decisive advantage 
of the variational method [IOllTT| is, that it extracts more 
information from the system by analyzing more than a 
single correlator. Again we use our example of the nu- 
cleon to illustrate this point. It is easy to check that the 
following two interpolators have the quantum numbers of 
the nucleon (C denotes the charge conjugation matrix), 

Oni{x) = eabc[ulix)Cj5db{x)]ucix) , (8) 

On2{x) = eabc[u^{x)Cdb{x)]j5Uc{x) . (9) 

Consequently both of them should give rise to Euclidean 
correlators that can be used to compute properties of the 
nucleon and of its excitations. 

We illustrate the effect of using different interpolators 
in Fig. [T] (taken from |8]), where we compare Euclidean 
correlators for different interpolators with the quantum 
numbers of the pion. The correlators are plotted as a func- 
tion of the Euclidean time t, normalized such that they are 
1 at t — 1. Since we use a logarithmic scale on the ver- 
tical axis, the exponential factor exp(— tM) gives rise to 
straight lines with slope —M. It is obvious, that beyond 
t — 5 the different correlators all show the same slope 
which corresponds to the mass M of the ground state. 
However, for small t the correlators have a different ad- 
mixture of exponential terms exp(— iM'), where M' > M 
is the mass of an excitation with the same quantum num- 
bers. If an interpolator couples strongly to such an exci- 
tation, this produces a steeper slope at small t. 

The central idea of the variational method is to use not 
only a single one of the possible correlators, but a whole 
r X r matrix of correlators, 



where each of the interpolators Oi, i — 1,2, ...r has the 
quantum numbers of the channel one is interested in. It 
is obvious, that such a correlation matrix extracts more 
information from the system than a single correlator. 

The determination of the physical observables from the 
correlation matrix can be cast into an elegant form |llj . 
One considers the generalized eigenvalue problem 



C(i)t;(") = A(i)(") C{to) «(") 



(11) 



where t^ <tis a, timeslice used for normalization. Explor- 
ing the particular form of the spectral representation of 
the correlation matrix (compare Eq. ([3])), 
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one can show, that the ordered eigenvalues A*^^^ > A'^^-' > 
A*^'^^ > ... behave as 

A(")(t) = e-(*-*«)^^" [l + e)(e-(*"*«)'^")] , (13) 

where M„ is the mass of the n-th state and Z\„ the dis- 
tance of Mn to the neighboring mass. Thus the ground- 
and excited states are disentangled, and each mass ap- 
pears in an individual eigenvalue. The largest eigenvalue 
corresponds to the ground state mass Mi, the second 
largest eigenvalue to the mass M2 of the first excitation 
et cetera. 

The fact that in Eq. \l'i\ the ground- and excited states 
are separated into individual channels allows for simple 
two-parameter fits of the eigenvalues. The corresponding 
ground- and excited state masses can be extracted in a 
stable and transparent way. 

We illustrate the power of the variational method in 
Fig. [21 where we plot the eigenvalues of the correlation 
matrix, again showing the example of the pion already 
used in Fig.[TJ The figure clearly demonstrates that, when 
plotted on a logarithmic scale, the slopes for the eigenval- 
ues are different, corresponding to the different masses of 
the ground- and excited states. 
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Fig. 2. Eigenvalues of the correlation matrix as a function of 
the Euclidean time. The data are again for the example of the 
pion, already used in Fig. 1. 



3.2 Constructing a basis of interpolators 



As any variational method, also the application to the 
excited states problem can only be as good as the basis of 
interpolators Oi one uses. Good interpolators should have 
several properties: 1) They should generate states from the 
vacuum that are as orthogonal to each other as possible; 2) 
These states should have a large overlap with the physical 
states; 3) The construction of the interpolators should be 
such that they can be implemented numerically in a cost- 
efRcient way. 

Several remarks are in order here: The true structure 
of the wave function of a physical state is unknown. Thus 
for constructing interpolators with a "large overlap" with 
the true physical wave function, only an educated guess, 
e.g., based on models is possible. However, even with a 
poor guess the method gives rise to the correct result, 
but the quality of the data might be poor. In practice, 
typically a large set of interpolators is implemented and 
subsequently the set of interpolators is reduced such that 
only the combination with the best signal to noise ratio 
is kept. If, e.g., an interpolator couples only very weakly 
to physical states, it will mainly contribute noise to the 
data, and should be left out in the analysis. An important 
consistency check of the method is to compare the results 
from different choices of the interpolators and to ensure 
that they agree within error bars. 

What are the different building blocks we can use for 
constructing hadron interpolators? The example of the nu- 
cleon, given in Eqs. dH), (IH), illustrates that different Dirac 
structures may be used to obtain different interpolators 
with the same quantum numbers. This freedom typically 
gives rise to two or three (depending on the channel) dif- 
ferent hadron interpolators with the same quantum num- 
bers /, J, P and C (where applicable). The different Dirac 
structures are an important ingredient, but obviously not 
enough for analyzing the tower of all excitations. It is gen- 
erally believed, that for a proper description of excited 
hadrons non-trivial spatial wave functions are essential. 



Recently the lattice community started to systemat- 
ically explore the possibility of implementing non-trivial 
spatial wave functions [3J-0 in hadron interpolators. Two 
different approaches for constructing such wavefunctions 
on the lattice are used: 

Basak et al have constructed large sets of possible 
hadron interpolators by using quark fields on lattice sites 
displaced relative to each other, typically involving near- 
est or next-to- nearest neighbors [7]. The operators were 
classified with respect to irreducible representations of the 
symmetry group of the hyper-cubic lattice, and in in this 
way the quantum numbers were assigned. 

A different strategy was followed in the work of the 
Graz-Regensburg group [3]-[5], [5]. This approach is based 
on so-called Jacobi smeared quark sources [91 , which is es- 
sentially a gauge covariant diffusion process, that leads to 
an s-wave type of wave function for an individual quark. 
Combining Jacobi-smeared sources of different width al- 
lows for radial wave functions with nodes. Application of 
an additional covariant derivative gives rise to p-wave type 
wave functions [HI . 

It is important to understand that the different inter- 
polators one uses in the variational approach are just the 
basis one offers to the system. The relative weights of the 
basis elements are not prescribed, but come out of the vari- 
ational calculation. Once the generalized eigenvalue prob- 
lem (|lip is solved and the eigenvectors v'"', n = 1, 2 ... r 
are known, one can define new interpolators O^^' as linear 
combinations of the original interpolators Oi , 



0(n) = J2 
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with the complex conjugate (denoted by *) entries v^ 
of the n-th vector as coefficients. The new interpolators 
O*^"-* are optimal within the given basis of interpolators Oi 
in the sense that they give rise to orthogonal correlation 
functions, 



(O(™)(t)o(«)t(o)) = A"(i)^„,. 



(15) 



as can be shown in a few lines of algebra. Thus the varia- 
tional method determines through the eigenvectors which 
linear combinations of the basis interpolators best describe 
a physical state. In principle the eigenvectors can be func- 
tions of the Euclidean time t. It was observed, that the 
eigenvectors are essentially independent of t in the range 
where the eigenvalues are dominated by a single exponen- 
tial. In Fig.[3]we illustrate this behavior for the eigenvector 
corresponding to the ground state, again using the pion 
example. Comparison with Fig. 2, where the correspond- 
ing eigenvalues are plotted, shows that indeed the ground 
state eigenvalue falls on a straight line up to i = 12, ex- 
actly the range where the eigenvector shown in Fig. 3 dis- 
plays plateaus for its entries. The property that also the 
eigenvectors show a signal stable in t is an important tech- 
nical tool for identifying the different states and for the 
determination of fit ranges. 
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Fig. 3. Entries of the first eigenvector of the generahzed eigen- 
value problem for the pion example already used in Figs. 1,2. 
The entries are plotted as a function of the Euclidean time t. 




Fig. 4. Final results for excited mesons taken from [5]. The 
horizontal bars are the experimental data, the symbols corre- 
spond to the lattice results. All masses are given in units of 
the Sommer scale ro = O.bfm. 



4 Results and future challenges 

4.1 Some selected results 

Having outlined the ingredients of an excited state calcu- 
lation in lattice QCD, let us now present some selected re- 
sults. We begin with discussing the outcome of a quenched 
calculation for excited mesons [5 . The calculation is based 
on meson interpolators constructed with the Jacobi smear- 
ing techniques outlined in the last section. The quark 
propagators were computed for several different values of 
the light quark masses and the results for the meson spec- 
tra were extrapolated to the chiral limit. The scale was set 
with the Sommer parameter and the value of the strange 
quark mass was determined from the K~^ meson. 

In Fig.Uwe show the results for the mass spectrum in 
different meson channels. The bars are the experimental 
numbers and the symbols are the lattice results. The errors 
for the lattice data are statistical errors. The plot shows 
that for pseudoscalar and vector mesons the lattice results 



agree well with the experimental values. The ground state 
masses typically are correct within 5 percent (10% for the 
/9-meson) and one or two excited states can be calculated 
which for the first excitations agree with the experimen- 
tal values within error bars. For the axialvector/tensor 
mesons the results are less convincing and we believe that 
here the set of basis interpolators is still not rich enough. 
Finally for the scalar channel the lattice data seem to co- 
incide with the experimental mass of the first excitation. 
Understanding the true nature of the scalar meson ground 
state is an interesting story of its own ~ it could, e.g., be a 
tetraquark state - and for a snapshot of the ongoing lat- 
tice work on this problem we refer the reader to the recent 
review |13) . 

Let us now come to discussing some results for baryons 
[5]. This quenched calculation is again based on the vari- 
ational method with basis interpolators constructed with 
Jacobi smearing. As for the meson case, the light quark 
masses were extrapolated to the chiral limit and the strange 
quark mass was set using the K~^ as input. 

Fig. [5] shows the results for positive parity in the Ihs. 
plot, while the rhs. is for negative parity. For both parities 
the ground state masses come out reasonably well - typi- 
cally with less than 5 percent deviation from the experi- 
mental numbers. As for the meson case one is wondering 
at this point why the quenched approximation is working 
so well. An exception is the positive parity A^'^ ground 
state which typically comes out about 20% too high [2]. 
A completely different story are the excited positive par- 
ity states. They come out typically 25% too high. Several 
possible reasons for this failure to describe these Roper- 
like states have been put forward, such as the quenched 
approximation, too small spatial lattice volumes and the 
fact that the light quark masses used still might be to 
large to capture the chiral dynamics necessary to get the 
masses right [2 . In general the opinion is, that the Roper 
puzzle will be solved only when dynamical simulations 
with light quarks on large volumes are completed. This 
is certainly an expensive enterprise, but the whole lattice 
community is pushing in this direction, and once these 
data are available, they will certainly be used for study- 
ing excited hadrons. 

To conclude the discussion of the baryons with a more 
positive aspect, we stress that the negative parity states 
come out very nicely. The lattice typically produces results 
that are on top of the experimental numbers. For a better 
resolution of the splitting between the ground- and the 
first excited state the error bars will have to be reduced. 
Nevertheless the results are sufficiently reliable, such that 
the lattice can make a prediction for states in the S and 
i? channels, where the corresponding masses are not yet 
well established experimentally (filled symbols in Fig.[S]). 



4.2 Upcoming challenges 

Let us finally come to the challenges that are waiting for 
excited state calculations on the lattice in the near future. 
We have already stressed the point that the quenched ap- 
proximation is the big unknown in the current results for 
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Fig. 5. Final results for excited baryons taken from [6]. The horizontal bars are the experimental data, the symbols correspond 
to the lattice results. In the Ihs. plot we show the results for positive parity baryons, while the rhs. is for negative parity. 



excited hadrons. We have also addressed the problem with 
the mixing to scattering states which has to be faced when 
light dynamical fermions are considered. In principle the 
techniques for disentangling scattering and bound states 
are ready to go [Mj : For scattering states the relative mo- 
mentum is quantized in multiples of the Matsubara fre- 
quency 27r/L, where L is the spatial extent of the lattice. 
Thus the energies of scattering states show a power-like 
dependence on 1/L which is absent for bound states, and 
comparing the results on different volumes allows to iden- 
tify scattering states. However, this is an expensive strat- 
egy which so far has been tested mainly in low-dimensional 
models and scalar field theories [TTlllSj. 

Another interesting challenge is the evaluation of ma- 
trix elements for excited states. Here some conceptual 
problems still have to be addressed, but in principle it 
should be possible to formulate this problem within the 
variational approach. A first appetizer of what the lattice 
could achieve in this direction is the calculation of the 
decay constant for 7r(1300) presented in [16] . 
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